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THE  IMPACT  OF  THE  THREE-WAVE  INSTABILITY  ON  THE 
SPIRAL  LINE  INDUCTION  ACCELERATOR 


I.  INTRODUCTION 

It  has  long  been  known  that  the  strong  focusing  magnetic  fields  that 

are  used  to  transport  the  high-current  beam  in  the  spiral  line  induction 

accelerator  (SLIA)^,  may  admit  instabilities  such  as  the  three-wave 
2-4 

instability.  The  SLIA  includes  solenoidal  focusing  in  the  straight 

sections  and  both  solenoidal  and  stellarator  (helical  quadrupole)  focusing 

in  the  curved  sections.  The  stellarator  fields  are  included  to  focus  the 

space  charge  of  the  beam  over  a  large  energy  bandwidth. 

Recent  results  have  shown  that  there  are  several  stability  regimes  for 

3  A 

the  three-wave  instability  ’  and  that  the  SLIA  proof-of-concept  experiment 

(PoCE)  can  operate  entirely  in  the  stable  regime.  What  has  been  less  clear 

is  the  significance  of  the  instability  as  the  beam  energy  approaches  oO  MeV. 

In  the  present  paper  we  will  address  the  difficulties  associated  with 

designing  a  three-wave  stable  SLIA  with  a  peak  beam  energy  of  50  MeV. 

The  three-wave  instability  is  a  resonant  interaction  between  the  beam 

centroid  motion  in  the  external  magnetic  fields  and  the  TE^  mode  of  the 

beam  pipe.  The  external  fields  include  the  axial  field  B  =  B  and  the 
r  v  z  o 

quadrupole  field: 

B  =  -B  k  [x  sin(k  z)  -  y  cos(k  z)]  (la) 

qx  q  q  q  q 

B  =  B  k  lx  cos(k  z)  +  y  sin(k  z)J  ,  (lb) 

qy  q  q  q  * 

where  B  k  is  the  quadrupole  gradient  and  k  is  the  quadrupole  wavenumber, 
q  q  q 

For  y  >>  1  the  stability  condition  of  Ref.  4  may  be  written  as 


K  > 
o 


K 


2um2 


4Yk2K  2  ,1/2 

+ _  q  q - 1 

yil(kq  -  U11/Y)J 
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(2) 


1 


2  -1/2 

where  y  =  (1  -  3  )  is  the  usual  relativistic  factor  for  the  beam, 

2  2 
K  =  eB  /y|3mc  is  the  cyclotron  wavenumber,  K  =  eB  /y(3mc  ,  w...  =  1.84/r  is 
oo  9  9  n  8 

the  cut-off  frequency  of  the  TE^  waveguide  mode  and  r^  is  the  radius  of  the 

1/2 

beam  pipe.  Note  that  for  <<  (k^u^j/4y)  ,  Eq.  (2)  reduces  to  the 

stability  condition  of  Ref  3.  For  example,  with  y  =  90,  =  0.048  cm~^  and 

r  =3  cm,  parameters  typical  of  a  high  energy  SLIA  bend^,  the  stability 

o 

condition  of  Ref.  3  is  valid  only  for  very  small  quadrupole  gradients, 
B  k  «  67  G/cm. 

q  q 

The  stability  condition  Eq.  (2)  gives  a  lower  limit  on  Kq  for  given 

values  of  y,  k  ,  B  k  and  r  such  that  as  the  beam  energy  increases,  greater 

values  of  Bq  >  0  are  required  for  stability.  With  the  complete  set  of 

stability  conditions  of  Ref.  4,  we  obtain  a  stability  diagram  as  shown  in 

Fig.  1.  Here,  stable  and  unstable  regions  of  (BQ,k^)  parameter  space  are 

delineated  for  given  values  of  y,  B  k  ,  and  r  .  In  Fig.  1,  y  =  7, 

Q  o 

B  k  =  200  G/cm  and  r  =  3  cm.  Note  first  that  the  SLIA  PoCE  will  operate 

q  q  g 

in  the  stable  region  that  appears  for  Bq  >  0  in  Fig.  1.  Note  also  that  in 

our  convention,  k^  >  0,  indicating  a  stellarator  field  with  a  right-handed 

pitch.  For  a  left-handed  stellarator  field,  k^  <  0,  stable  operation 

requires  B  <0. 
o 

The  difficulty  of  satisfying  Eq.  (2)  at  high  energy  can  be  seen  in  the 
stability  diagram,  shown  as  Fig.  2,  for  y  =  90,  B  k  =  73.4  G/cm,  and  r  =  3 

MM  © 

cm.  Figure  2  shows  that  Bq  >  8.5  kG  would  be  required  for  stability.  This 

is  confirmed  by  numerical  solutions  of  the  dispersion  relation.  Because 

such  high  values  of  Bq  may  not  be  desirable  and  because  three-wave  growth 
-1/2 

rates  fall  as  y  ,  it  may  be  favorable  to  operate  in  a  mildly  unstable 


regime  at  high  energy.  Typical  mildly  unstable  parameters  for  a  y  =  90  SLIA 
bend  are  shown  in  Table  I.  Our  uncertainty  as  to  how  much  growth  is 


2 


tolerable  suggests  to  us  that  we  consider  alternate  strategies  for  reducing 
or  eliminating  growth  in  this  regime. 

In  the  present  paper  we  discuss  three  such  strategies.  In  section  II 
we  discuss  detuning  the  instability  by  varying  the  axial  guide  field. 
Section  II  contains  a  discussion  of  the  possibility  of  damping  the 
instability  through  the  use  of  a  low-Q  beam  pipe.  In  section  IV,  we 
consider  the  possible  use  of  a  discrete  alternating  gradient  quadrupole 

focussing  system,  also  known  as  a  FODO  lattice.  This  opens  up  the 

possibility  of  detuning  the  instability  by  varying  the  quadrupole 
wavenumber.  We  summarize  our  results  in  section  V. 

II.  DETUNING  IN  THREE-WAVE  UNSTABLE  REGIME  II 

Recent  ELBA**  simulation  results  showed  saturation  of  the  instability  in 

three-wave  unstable  regime  II  (k  /2  +  2K  <  K  <  K  ).  Unstable  regime  II 

q  q  o  c 

corresponds  to  unstable  region  II  of  (BQ,k^)  parameter  space  as  shown  in 

Figs.  1  and  2.  This  saturation,  which  occurred  without  emittance  growth  or 

current  loss,  was  the  result  of  detuning.7  In  this  scenario,  the  beam  loses 

energy  to  the  TE^  waveguide  mode  and,  at  the  same  time,  the  parallel 

velocity  of  the  beam  is  reduced  as  the  transverse  motion  increases.  Both  of 

these  effects  increase  the  cyclotron  wavenumber  K  such  that  the  three-wave 

o 

interaction  is  detuned.  This  detuning  changes  the  frequency  <*>  of  the 
unstable  mode.  The  transverse  beam  centroid  motion,  which  has  wavenumber 
w/Pc,  cannot  respond  promptly  to  this  frequency  shift  with  the  result  that 
the  beam  motion  is  damped. 

Diagnostics  of  beam  energy  and  axial  and  parallel  velocity  of  the  beam 
centroid  from  several  ELBA  runs  that  exhibited  this  saturation  showed  that 
both  y  and  0z  decreased  by  several  percent  as  the  instability  grew.7  This 
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shifted  the  cyclotron  wavenumber  by  =10%  from  the  point  when  the  instability 
began  to  grow,  which  was  taken  to  be  one-third  of  the  saturation  distance, 
to  the  point  of  saturation.  Solutions  of  the  dispersion  relation  for  those 
parameters  showed  that  for  a  10%  change  in  either  Bq  or  y  such  that  Kq 
increased  by  10%,  the  unstable  interaction  point  in  («,k)  space  shifted  so 
as  to  stabilize  the  previous  region  (in  frequency  space)  of  peak  growth. 

Analysis  confirms  that  detuning  in  Kq  is  effective  in  three-wave 
unstable  regime  II  but  has  little  or  no  effect  in  regime  I  (Kc  <  k^/2  -  K^). 
To  see  this,  we  must  reexamine  the  three-wave  interaction.  Here  we  consider 
unstable  waveguide  modes  with  w  >  0,  but  the  discussion  carries  easily  to 

the  w  <  0  case.  As  one  sweeps  through  the  possible  values  of  Kq,  the 
equation  of  the  unstable  beam  mode  in  (w,k)  space  in  the  1^  ■*  0,  K  ■*  0 

4 

limit  changes.  For  Kq  <  0  (in  unstable  regime  I),  the  equation  of  the  beam 
mode  is 

k  =  «/c  -  kq  .  (3a) 

For  0  <  Kq  <  k^/2  (also  in  regime  I),  the  unstable  mode  is 

k  =  w/c  +  K  -  k  .  (3b) 

o  q 

For  kq/2  <  Kq  <  kq  (unstable  regime  II),  the  unstable  mode  is  again  given  by 

k  =  »/c  +  K  -  k  .  (3c) 

o  q 

For  K  >  k  ,  the  beam  is  stable  in  the  I,  -»  0,  K  -»  0  limit, 
o  q  b  q 

We  see  that  in  unstable  regime  I,  shifts  in  Ko  have  little  effect 
because  the  beam  mode  is  is  either  independent  of  Kq  or  restricted  to 
Kq  <  kq/2.  In  unstable  regime  II,  detuning  is  possible  and  is  observed. 
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These  results  suggest  that  it  may  be  possible  to  operate  with  low 

growth  rates  in  three-wave  unstable  regime  II  by  employing  cyclotron 

detuning.  Here  one  periodically  shifts  the  magnitude  of  the  axial  guide 

field  to  shift  the  instability  from  frequency  to  frequency,  thus  keeping 

transverse  beam  motion  at  any  given  frequency  below  significant  levels. 

Numerical  solutions  of  the  dispersion  relation  show  that  the  latter  two 

cases  of  Table  I  can  be  detuned  with  shifts  in  B  of  <10X. 

o 

A  refinement  of  this  detuning  technique  might  include  tuned  stubs  at 
the  location  of  each  change  in  the  axial  magnetic  field  value.  In  theory, 
these  could  be  adjusted  so  as  to  reflect  the  low  level  TE^  mode  at  the 
point  where  the  unstable  frequency  is  shifted  and  a  new  mode  begins  to  grow. 

Note  that  Eqs.  (3a-c)  show  that  detuning  the  instability  via  changes  in 
kq,  rather  than  Ko,  may  be  possible.  Results  for  alternating  gradient 
quadrupole  focusing  systems,  presented  in  section  IV  below,  confirm  this 
possibility.  Fabricating  a  helical  quadrupole  magnet  with  a  varying  pitch 
length,  however,  is  considered  prohibitively  expensive. 

III.  DAMPING  VIA  LOV  Q 

Because  growth  rates  in  the  high  energy  regime  are  small,  it  may  be 
possible  to  damp  the  instability  by  lowering  the  Q  of  the  beam  pipe.  If  we 
assume  that  the  presence  of  damping  does  not  change  the  nature  of  the 
instability,  the  three-wave  interaction  will  be  stabilized  when 

Q  <  w/2 T  (4) 

*—  OD 
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where  «  is  the  frequency  of  the  unstable  waveguide  mode  and  is  the  growth 

O 

rate  from  the  linear  dispersion  relation  with  Q  =  ®.  Analytical  results  in 
which  damping  was  added  to  the  dispersion  relation  have  verified  Eq.  (4). 

For  example,  for  the  three  parameter  sets  of  Table  I,  the  instability 
can  be  damped  with  the  rather  high  Q  values  of  850,  3700  and  40,000 
respectively. 

Another  example  is  the  NRL  modified  betatron  accelerator  (MBA)  in  which 
the  addition  of  strong  focusing  fields  has  aided  the  successful  acceleration 
of  a  1  kA  electron  beam  through  >38000  turns  in  the  100  cm  major  radius 


device 

to 

obtain  15  MeV.^ 

Typical  MBA 

parameters  are 

beam  current 

i.  =  1 

kA, 

B  =  -2  kG,  B  k  = 

20  G/cm,  k  = 

0.06  cirT^  and  r 

=  15  cm.  As 

b 

o  q  q 

q 

g 

the  beam  is  accelerated  from  5  to  15  MeV,  the  growth  rate  varies  from 
TJc  =  5.10  x  10  ^  cm  *  to  3.04  x  10  ^  cnT^.  Because  the  frequency  of  the 
unstable  mode  is  insensitive  to  changes  in  the  cyclotron  wavenumber  in  this 
regime,  the  unstable  frequency  remains  constant  at  t*>/c  =  0.16  cm-1  (w  =  4.8 
GHz).  Here  Q  <  150  is  required  for  stability.  The  MBA  vacuum  chamber  has 
Q  =  150.  Three-wave  stability  was  confirmed  by  radiation  measurements  using 
a  broad-band  antenna  which  showed  radiation  levels  not  exceeding  those 
expected  from  single-particle  emission. 

IV.  DISCREET  QUADRUPOLE  FOCUSING 

Discrete  quadrupole  focusing  systems,  also  •'ailed  F0D0  lattices,  have 
stability  properties  that  differ  from  helical  quadrupole  systems. ^  The 
focusing  field  for  such  a  system  may  be  approximated  by 

B  =  -B  k  cos(k  z)y  ,  B  =  -B  k  cos(k  z)x  (5) 

qx  q  q  q  qy  q  q  q 
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where  k  is  the  FODO  wavenumber,  B  k  is  the  FODO  gradient,  z  is  the  axial 
q  q  q  6 

coordinate  and  x  and  y  are  transverse  coordinates.  Note  that  the  FODO 


lattice  field  has  no  preferred  direction  so  that  we  expect  Bq  <  0  to  be 

equivalent  to  Bq  >  0.  In  fact,  when  the  dispersion  relation  is  solved  and 

approximate  stability  boundaries  are  found**,  we  obtain  the  stability 

diagram  shown  in  Fig.  3.  Here,  y  =  7,  B^k^  =  200  G/cm  and  r^  =  3  cm. 

It  is  clear  from  the  stability  diagram  that  the  overall  stability 

properties  are  not  as  favorable  as  in  the  helical  quadrupole  case.  However, 

we  find  that  for  identical  values  of  I, ,  y,  B  k  and  k  ,  peak  growth  rates 

o  q  q  q 

1/2 

are  lower,  typically  by  2  ,  in  the  FODO  configuration.  Furthermore,  the 

flexibility  of  the  FODO  system  allows  growth  rates  to  be  reduced  via 
detuning  in  the  FODO  wavenumber  k^. 

We  find  that  the  FODO  wavenumber  can  be  changed  so  as  to  shift  the 
unstable  mode  from  frequency  to  frequency.  Also,  this  detuning  mechanism  is 
applicable  to  low  axial  magnetic  fields  (0  <  Kq  <  k^  -  2k^)  wbere  cyclotron 

detuning  is  not  effective.  This  is  demonstrated  in  Fig.  4,  which  shows 


growth  rate  versus  frequency  for  a  FODO  lattice  with  Bq  =  0,  B^K^  =  200  G/cm 

and  two  different  values  of  k  =  0.4  and  0.6  cm-*.  Here,  I.  =  10  kA,  y  =  7 

q  b 

and  r  =3  cm.  We  see  that  the  unstable  range  at  k  =  0.6  cm-*  is 

g  q 

completely  stabilized  when  the  FODO  wavenumber  is  shifted  to  k^  =  0.4  cm  *. 
As  with  the  helical  quadrupole  case,  detuning  will  be  effective  if  the  beam 


centroid  wavenumber,  w/3c,  is  changed  often  enough  that  significant 


transverse  beam  motion  does  not  develop.  This  scenario  may  be  complicated 
by  the  multiple  unstable  wavenumbers  associated  with  the  unstable  frequency 
of  the  FODO  system,  making  it  difficult  to  detune  the  TE^  mode  both  in 
frequency  and  wavenumber  as  can  be  done  in  the  helical  quadrupole  case.  We 
note,  however,  that  the  importance  of  this  is  not  well  understood. 
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Furthermore,  de  igning  such  a  minimum-growth  FODO  system  may  be  somewhat 
time  and  computing  intensive.  Note  that  the  instability  in  this 
configuration  has  been  verified  using  the  ELBA  simulation  code.** 


V.  CONCLUSIONS 

Previous  results  have  shown  that  three-wave  stability  in  the  SLIA  is 

easily  maintained  for  energies  up  to  ^0  MeV.  We  find  that  growth  rates  for 

-3  -1 

energies  exceeding  40  MeV  can  be  <  10  im  It  may  be  possible  to 

-3  -1 

determine  whether  growth  rates  of  order  10  cm  are  permissible  over  tne 
last  2-3  bends  of  r.  50  MeV  SLIA  via  PoCE  experiments.  If  such  growth  rates 
are  not  acceptable,  we  have  o  tlined  three  strategies  for  stabilizing  or 
reducing  the  three-wave  interaction.  These  are  1)  cyclotron  detuning, 
wherein  the  instability  is  detuned  over  the  length  of  a  given  bend  by 
varying  the  axial  magnetic  field,  2)  damping  the  instability  through  the  use 
of  a  "low-Q"  beam  pipe  and  3)  discrete  quadrupole  focusing,  for  which  it  is 
feasible  to  detune  the  instability  by  varying  the  discrete  quadrupole 
wavenumber  over  the  length  of  the  curve.  We  show  that  "low-Q"  can  include  Q 
as  high  as  40,000.  Note  that  de-Qing  and  cyclotron  detuning  are  most 
effective  in  thee-wave  unstable  regime  II,  where  the  interaction  frequency 
is  highest  and  is  sensitive  to  the  cyclotron  wavenumber. 
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Table  I 


Bq  (kG) 

B  k  (G/cm) 
9  q 

_i 

k  ( cm  * ) 

q 

co/c  (cm  *) 

T/c  (cm  *) 

-6.91 

235.0 

0.064 

3.41 

2.0  x 

10  3 

7.10 

73.4 

0.048 

11.23 

1.5  x 

io-J 

7.71 

58.3 

0.048 

28.56 

3.5  x 

10'4 

Unstable  growth 

rates  and 

frequencies  for 

typical  SLIA 

achromat 

parameters 

with  y  =  90,  1^ 

=  10  kA,  r 

g 

=  3  cm  and  bend 

radius  R  * 

0 

250  cm. 
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Fig.  1  —  Stability  diagram  for  beam  energy  7=7,  quadrupole  gradient  Bq  k,  =  200  G/cm  and  waveguide  radius  rg 
cm 
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Fig.  2  —  Stability  diagram  for  beam  energy  y  =  90,  quadrupole  gradient  Bq  =  73.4  G/ cm  and  waveguide  radius  r?  = 
3  cm 
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Fig.  3  —  Stability  diagram  for  a  FODO  lattice  with  beam  energy  7  =  7,  quadrupole  gradient  Bq  ks  =  200  G/cm  and 
waveguide  radius  rg  =  3  cm 


w/c  (cm  ') 

Fig.  4  —  Growth  rate  T  versus  frequency  ui  for  a  discrete  quadrupoie  focusing  system  with  B0  =  O.  Bq  =  200  G/cm, 
Ib  =  10  kA,  y  =  7,  rg  =  3  cm  and  two  different  values  of  k,,  =  0.6  cm'1  (solid)  and  0.4  cm'1  (dashed). 


